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ABSTRACT 


Geomagnetic micropulsations were continuously 
recorded during the summer of 1968 at the University of 
Alberta Geophysical Observatory. Simultaneous measure- 
ments were made at a series of field hearse in Alberta 
WLeh truckemoennvedc instrumentation py The sanalogue. records 
were electronically converted to digital form and power 
spectral density estimates were obtained using Fast Fourier 
Transform techniques. From these spectral estimates the 
polarization characteristics were computed for the horizon- 
tal magnetic field variations having a period range from 
MOmtee27s)ceconas. sihe polarization: characteristics are 
displayed as a function of period and time of day. 

The averaged results indicate that the magnetic 
Sienalewas DOlarizedeirom north vo NOO°E. Micropulsations 
falling into the Pie class were found to be predominantly 
polarized in the counterclockwise sense, while Pils were 
polarized in the clockwise-sense. Pe3's were polarized in 
the clockwise sense in the local morning and in the counter- 
clockwise sense in the local afternoon while Pc4's were 
counterclockwise before local noon and clockwise in the 


afternoon. 
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CHAPTER 1 


INTRODUCTION 


ea) Definition of Micropulsations 


For over a century it has been known that the 
Rarth's geomagnetic field exhibits amplitude variations. 
These variations can range in period from small fractions 
of a second to many millions of years. It has been 
realized that mueh information concerning the Earth's 
interior and crust, ionosphere and magnetosphere can be 
obtained from studies made on these variations. In recent 
years interest has been shown in the magnetic field 
fluctuations known as geomagnetic micropulsations. This 
class includes all the PPuctustlens having a period 
falling in tne range from aoouv 0.2 seconds. to 10 minutes. 
Their amplitudes can vary from the smallest detectable 
value (which is the order of a few milligammas at present) 


up to a few hundred gammas on rare occasions. 


Studies of the characteristics of micrepulsations 
have led to the separation of this activity into various 
groups. Many authors have devised their own classification 
ayatems, but.a.certain,amount.of difficulty is encountered 
due to uncertainties in the origin of micropulsations and 


their inter-relationships. To standardize the nomenclature 
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used, a classification system was approved at the thirteenth 
General Assembly: of the 1UGG in Berkeley, California 


ClAcolewecra Irs ods 


in this classification, micropulsations were 
G@ivided into, two main elasses: those of regular and 
continuous appearance ealledqvrc, and those Gen anearregyular 
fPormecaiicad sy... The: Pe group was further divided into, five 
SUMPSrOUpe sana tie Fiseroup into two subgroups. These 


SvOsrOUpSsaremsnowhs in. Table. 1. 


TABLE I 
CLASSIFICATION OF MICROPULSATIONS 


CLASSIFICATION PERIOD CLASSIFICATION PERIOD 
(SECONDS ) (SECONDS ) 


0.2-5.0 

5.0-10.0 

10-45 
eeu 


150-600 





A general picture of the pulsation spectrum is shown in 


Figure 1.1 (after Campbell, 1966). 
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ie ce hetey | Mal, General picture of the geomagnetic pulsation 
Spectrum for the lower frequencies (Campbell, 


1966). 
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Hee Literature Review 


Numerous studies have been made to gain further 
insight into the. source mechanisms of micropulsations. 
Considerable interest has been shown in recent years in 
their polarization characteristics. This thesis analyses 
these characteristics. for.micropulsations having, a-period 
Pens ee Ole Lm Onek MS eCONdS. «4 ounce nDOLNAPIpamdsFc 
pulsations were.analysed, it-is.necessary tocdetermine how 
pOlangization, characteristics vary.between-Pil-2's.+and 


Fo3-5i's. 


Polarization studies have often been done by 
tracingPoueethe-+tip-of’ thevmagnetic™"perturbation*vector in 
a plane, (the hodogram technique). If the variations of 
the magnetic Ffileld*®are polarized*thens"in?’general, an 
eiiiptical formewiltl*®be traeted*out?~+~Linéar"and*eireular 


Prormsreceursenly @rarely. 


Since polarization properties are heavily 
dependent. on the time of day, latitude, and frequencies 
present in the micropulsation event, a brief review of 
observations on these factors is presented here. Since 
complete reviews are to be found elsewhere (Jacobs and 
Westphal (1964), Campbell (1967), and Saito (1969)), it 


was felt necessary to include here only observations on the 
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Perlodpandstamesote daynoccurrence of.~micropulsatvions, 


te eral Occurrence Frequencies and Period Content 


ofalMicropulsations 
(a) Pi Micropulsations 


Pie pulsations have been found to predominate in 
the nighttime and have a maximum occurrence frequency 
around local midnight, (Romand and Cardts, 1962). 
However, they have also been detected in the daytime 
(Voelker, 1968). Pil pulsations often appear at the 
beginning of Pi2 events and also appear as riders on the 
longer period Pi2's. Jacobs and Sinno (1960a) have 
shown that their amplitude falls off to the north and 
south of the auroral zone. It has been found by 
Yanagihara (1959) that, as solar activity increases, Pi2 
occurrence decreases, but ites J UStetne Olpositen ror 
Pil pulsations. However, the Pid occurrence frequency 


increases with increasing magnetic activity. 


Pi2's can exhibit fairly complex frequency 
spectrums. Rostoker (1967a) found a linear relationship, 
in middle latitudes, between the number of frequency 


components in Pi2 events and Kp index. The number of 
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frequency components increases as the Kp index increases. 
On the average, he finds that the dominant periods lie 
between 100 seconds to 130 seconds. Rostoker also reports 
that these frequency components are not generally har- 
monics of a fundamental frequency so it is possible to 
consider that each frequency has its own separate gener- 
ation mechanism. However, Hirasawa and Nagata (1966) at 
the Kakioka Field Station (¢ = 36°14'N, A = 140°11'E) 
found’ that in’ 58 samples of Pi2%pulsations ,\ 74% consisted 
of a fundamental period and at least one harmonic 
component. Nwaigwe et al. (1967) analysed the pulsation 
activity recorded in south-west England in the period 
range from 10 seconds to 150 seconds. They reported an 
average period of about 80 seconds at night, but it 
varies from 110 seconds at 2230 GMT to about 40 seconds 
at O400 GMT. Saito and Matsushita (1968) found the 
period of Pi2 to be between 60 to 90 seconds. They 
obtained this result by using digital dynamic spectrums 
of pulsations recorded at Fredericsburg (6 = 38°12'N, 


No BLO 293364 WS . 


(b) Pe Micropulsations 


The Pe type pulsations have been found to be more 


prevalent in the daytime hours (Jacobs and Sinno, 1960b). 
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Pe2 and Pc3 pulsations usually have an amplitude of about 
On5 gammage Campbell be( 1967 psbelievées “Ghatcthe*duration of 
these signals is related to the sensitivity of the recording 
device. They are always present but are at times too small 
to detect. Pc2 and Pc3 are reviewed together here as there 
Seems to be no distinct spectrum separation between them. 
Hirasawa and Nagata (1966) report that Pc3 pulsations at 
the Kakioka Field Station appear only in the local daytime 
from 0600 hours to 2000 hours and are most active around 
local noon. Nwaigwe et al. (1967) found that during the 
mMonning thours’§ ePic3 wactiviity thad.tan “average “period yof about 
SS eeCCnds: andwintithe afternoon f#tanraverage “perFod of about 
40 seconds. This trend was also noted by Voelker (1968) in 
hisstudies of stations in Europe.. Also Hirasawa and 
Nagata (1966) noted that the average period increased from 
25) seconds tat 0600 hours “LT rto133 ‘seconds ‘atvwlL800 hours LT. 
These average periods are smaller than those found by 
Nwaigwe et al. (1967). This is perhaps explained by the 
fact, thaw Voelker (1960), finds 4. svystemauiconermeas cin thne 
period. af Pe? and Pe 3 sin whe wayuime, of sthe morthssouch 
component of the magnetic field with increasing latitude. 
Kato (G00) 2 on, thevotherm hand, staves that bcs souleations 
have the longest period around noon and the shortest period 


around dawn. Entirely different results were found by 
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Christoffel and Linford (1966) in New Zealand. They found 
the average period to decrease from 40 seconds at O400 
hours LT to 27 seconds near local noon, and then increase 
again to 40 seconds at about 1800 hours LT. These results 
were obtained by visual analysis. Some of this confusion 
Maye-ac lees GUevOMthel factitthats Pc3 pulsations with a 
shorter period at noon (U-Type) tend to appear in the 
Sunspot minimum years while those with longer periods 
around noon (inverted U-type) appear at the times of sun- 
Spot maximum. Pim aienra temper ou ot, bc sisunas been 

Found %0 increase with decreasing sunspot activity while 
Pc4 period decreases. Pc4 pulsations tend to be of the 
U-type and are active during sunspot minimum years. 
Therefore, Pc3's can be distinguished from Pc4's during 
Sunspot. maximum years, but during sunspot minimum years, 
Pes euememnacskedsoy Pe. 's Gntruding into. them. The 
CJUbneIeoer Lod svarlatlon observed will depend, then, on 


the year of the observation (Saito, 1969). 


Pc4 micropulsations have an amplitude that ranges 
from 5. Eow2O sgammas win stie higher Latitvudescwandsan 
eventsconulact® fromn,lO minutes to Sseyeral hours, They 
eanwalso take a beating form, or have a .dampedsappearance. 


Hirasawa and Nagata (1966) report that Pc4 pulsations are 
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COmvLinveus throughout the days but are most active in the 
morning hours with the maximum occurring around 0800 hours 
LT. They found that the average period was 50 seconds 
around O7OU ours: Ll ana lll seconds around 1900 hours LT. 
Saito (1964) found an occurrence maximum at midday in 


middle and equatorial latitudes. 


Pe5 pulsations have the largest amplitudes and are 
someulmes Called Pao"s or giant pulsations. They have 
quite a Sinusoidal appearance and can be damped, but if 
DicaitecweuMomecul vi ace hioi, “some distorc1on can occur. 
Pe>'s display a daytime occurrence pattern. Campbell 
(1967) states that near the northern auroral zone, Pc5 
pulsations occur most frequently near 0600 hours and 1800 
VOUtoM ee iOWCVeT.= ap middle latitudes, a noon maximum is 
also reported. Herron (1967) averaged 80 records from the 
Lebanon state, Forest ,,.NeweJersey,estation tos obtain the 
periods oOf- spectralopeaks:. tele) foundsthatr the] spectrum 
peaks at 227, 285 and 393 seconds. These spectral peaks 


had.-the highestyamplitudeion therdays ofg@highnkptindex? 


Diice one ee irVves=a Senerals review Of © UiLe 
characteristics of micropulsations. While many interesting 
facts have been obtained, much work still remains to be 


done. Many workers have analysed records that display 
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only simple | funcontaminated: signalw® They ittlave (avoided 
records that have several frequency components present and 
haVvewerre mula vamplie tudes, ~Thiswisvespectial iy true fon 


Ssudike aso fiepolarization. 


dec2 wre Polarization Characteristics of Micropulsations 


Polarization characteristics have been determined 
mainly by visual inspection of magnetograms and by plotting 
tiie enusepotme or Chesperturbation vector. However, Paulson, 
Egeland and Eleman (1964), and Paulson (1968) have reported 
Workmsdone by analyolcal techniques. A review is given here 


OMmeCte results EOUNG by Various researchers. 


(fae) Sense: or Polarization 


Zyoin (1907) reports that results obtained in the 
Wsskechowea that the sense Of rotation ef the polarization 
ellipse was 80 to 90% counterclockwise for all micro- 
pulsations in the period range from 30 seconds to 1200 
seconds in the local morning hours in middle northern lati- 
tudes. In the evening hours, however, 60 to 70% of the 
micropulsations had a clockwise sense of rotation. He 
reports that separation of the two senses occurs one or two 


hours before local midday - midnight. 
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Zyban. aid notwappear tosmake-a@ distinetion; between 
Pieenc) ve events. | Rosteker (1967b); however, studied, 20 
separate Pi2 events recorded simultaneously on a station 
array extending from Victoria to Montreal, Canada. He 
found that Pi2's are polarized predominantly in the 
counterclockwise sense atemicdle latitudes: in fhe morthern 
hemisonere.. he found that, even ifthe source of the Pic 
micropulsations lay in the center of his east-west station 
array, the sense of polarization almost always remained 
Ponecant, P kato et al... (1956) found, however, that Pi2's 
were polarized in the counterclockwise sense after local 


Micnieiieeene if che clockwise sense before, 


Onristorrel end Gintord (1966) réport that in 
New ezealaned 61% of =the Pi2 events: had a clockwise sense 
and 37% were linear. Here, one must take into account the 
predicted reverse sense of rotation between the two 


hemispheres. 


It seems, then that the preferred sense of rotation 
for Pi2's is still somewhat confused. The average results 
of all the local nighttime records analysed in this thesis 
showed that Pi2's are predominantly polarized in the 


counterclockwise sense. The Pil band however shows 
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airPrerent results as will be discussed’ in-Chapter III. 


in New Zealand, Christoffel and Linford (1966) 
found that 56% of Pc3's analysed between 0700 to 1900 
NOuUrse.) were polarized in the counterclockwise sense and 
11% were clockwise; the rest. were “Wineer., Mey reported 
that nighttime Pc3 and Pc4 events were weak and not a 
sufficient amount of data was collected to give signif- 
icant results. Campbell (1967) reported that Kawamura et 
Sveaello) sound that only 5% of Pe2—3 events had a linear 
form.) tne resy were elliptical. “These results were found 
in middle- latitudes. About 75%, of. the ellipses hada 


fet oworemiunor tO Mayor axis et Less than-0.>. 


Por sPe wis. Camppell (1967) states that the sense 
of polarization is predominantly counterclockwise in the 
northern hemisphere. Mather, Gauss et al. (1964) reported, 
however, that the sense of Pc4's was divided into clockwise 
during the day (0600-1800 hours LT) and counterclockwise 


ab night. (1300-0600 hours LT) in middle northern latitudes. 


One hundred and ninety Pc5 events were recorded at 
Point Barrow (¢ = 71°18'N, > = 156°46'W) ,College (@ = 64° 
be'W, A = LA7°50'W) and Sitka (4 = By CON Ne Me =i 5 2G WD 


by Kato and Utsumi (1964). The analyses indicated a 
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counterclockwise sense of rotation before local noon and a 
clockwise-.sense in the. afternoon. In the middle latitudes, 
Sanos (9680+ and? Saitoe(1964)),.report.thataPc5rpolarization 
seems to have four sectors in which the polarization 
changes alternately from clockwise to counterclockwise. 

The results obtained in this thesis suggest that events 
Clagsitied as Ped5's are polarized: in the counterclockwise 
sense. The longest period sampled, however, was 273 
seconds and these results are not really representative of 


the wntire Pc5 band. 


(b) Anele of Polarization 


ihe saneile, Ol pOlarigacion is defined ie tle ules. 
ao che angle thar the major|axis of DelerizarLom en lipse 
makes with geographic north. It is Gelgned cs POstu1ve =o 
it is to the east of north, and negative if it is to the 
west. Very few conclusive results seem to have been found 
previously on. this aspect of micropulsations. JZybin (1967) 
states that the angle depends mainly on the characteristics 
of each recording station and not on the type Grim cropuld— 
sation being received. He observes that Stations, near 
simple but large extended geological structures have Cie 


Maori. s.of the pobarization ellipse perpendicular to the 
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Surikesol the soructure., Zybin made, this statement with 
reference to hourly averaged results as individual 
measurements of the angle show a good deal of scatter. 
PhiSescauter, tnen,reflects. the, variations a Ee 
micropulsations themselves, whereas mean values. show the 


influence. of geologic structure. 


Table II, reproduced from Zybin (1967), shows the 


preferred angles of polarization at several stations. 


pA Biles PLD 


PREFERRED POLARIZATION DIRECTIONS OBSERVED AT 
DIFFERENT STATIONS 


AZIMUTH OF 
STATION PREFERRED RESEARCHERS 
DIRECTION 


Aburatsubo Hatakeyama (1938) 
Bermudas Santirocco and Parker (1963) 


Borok Kalashnikov and Zybin (1960) 
Zybin (1966) 


Gottingen Untiedt (ai 961) 
Kiruna Paulson et 41... (1965) 


Lovozero Barsoukov and Zybin (1961) 


Petropavlovsk Zybin (1966) 





Victoria DUrfus and’ Shand? (1958) 
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Zybin found that those stations located near regional 
extended geoelectric structures and shorelines 

(Aburatsubo (Japan), Borok (% = 58°02'*N, A = 38°58'E), 

GU@e eet Ceae od sa A = 9°50 hey Kiruna (o)= 67°50'h., 
Mee e0 25) ), and, Lovezero (¢ = 67°58'N, A = 35°58'E)) have 


eevke,crred, C1recuLon, OF pOlarizatiom perpendicular. te the 


SeiakKeow slandsstavions such as the. Bermudas, Petropavilovsk 
Cuero CON. er 150: SOLE), ands Vicroria (bo =" 40° 30'N, 
XW = 123°24'W) are much more difficult to interpret. Since 


the Rocky Mountain Range was within 200 miles of all the 
recor.dlipeouata Ons Uuscdslnethiis, toesis, comparisons were 
made between the angle of polarization and the strike of 


these mountains. 


iyo GlOG4)seeporreduthat., Stapistically, the major 
exicwOrpeenem polarization ellipse rotates, cleockwiservirom the 
everagce prevailing direction before midnight and rotates 
eounterclockwise in the hours after midnight for both Pi and 
Pe events. Around midnight, however, he found that the angles 
ace,close to their preferred direction.» From st tudes on 
mecels Zybing(l96/) found that the diurnal variation or tae 
angle is due to the change in the sense of polarization. 
Correlations were made in this thesis between these two 


Polarization aspects. 
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ZO ge ClOod Wal ceoer ounce pat befsshad. an. angle 
OpepOlearizallem theausrotates counterclockwise relative to 
Une average prevailing direction, in the,morning, while in 
ules evenine «it rotates, clockwise. Campbell (1967) states 
Unhepecthe polarization angle.is zero (1.6... the major axis 
of the ellipse is north-south) for Pe2's and Pe3's. He 
feels, however, that this is somewhat dependent on the 
Roc Viacine Ol day. a kato (1964) supports this and finds 
that the north-south component of the magnetic field 
disturbance vector predominates at low latitudes. This 
Mmeancacres nearly linear polarization, showing that the 
Prieto! (2x15 Or polarization is in the north-south 
Orcecrron. Lhe. .dturnal vartavwon ofthe major. axis: in 
Japan was found to be to the NNW during the day and to the 
NNE ae Mient., (Kate and Saito, 1959; Kawamura, Kurusu eét 
eel oo). Christoriel and Lindferd- (1906), on the other 
heandwd.e not seem to find any dominant angle for Pc3s'’s in 
New Zealand. For a given Pc3 event, the direction tended 
Toscemaghtairly constant in, the altennoon, butwea cradual 


change began near sunset. 


According to Campbell (1967), PceS pulsations are 
Uslialivy elliovtcally polarized. ~The angle of poleri7acien 


usually changes between the morning and evening peaks of 
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the diurnal frequency pattern (Kato, 1964). 


ike s Purpose ofimhes Research 


Many theoretical attempts have been made to explain 
the origin of micropulsations and consequently, have 
predicted a variety of polarization Characteristics. orlhe 
results presented in the previous section indicate that an 
entirely coherent experimental discription of the polariz- 
aulon characteristics has yet. vo be determined. “TA. large 
amount of statistical information about geomagnetic 


micropulsations must be accumulated. 


Until the last few years, polarization information 
was obtained by utilizing the hodogram technique. This 
technique réquires that- the signal be nearly per lod ler, 

Any contamination due to random noise or the superposition 
Of several signals at different periods will result=in an 

extremely confused hodogram. This was clearly demonstrated 
EPyepoauleon el 960). ) Thus, detailed polarization studies of 
complicated magnetic field variations have-not, in general, 


been attempted. 


Isis “the plrposesof this thesis vo present the 


results using power spectral analysis and. the techniques 
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outlined by Fowler et al. (1967) to analyse complicated 
magnetic: signeds yY'Tfhis+technique®was,appiteds’for the 
Hinscmorme, 60,a large-number of v~records to determine 
statistically the hourly averaged values of the polarization 
echanacteristics-as.a function, of periods and.bime,of day. 
huge polarization SS ere: of the geomagnetic 
micropulsation spectrum are clearly displayed for 


imcperpretati on.purposes.. 
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CHAPTER®2 


INSTRUMENTATION AND DATA PROCESSING 


aes 1. Signal Detection 


Vanlavionsi®in=theesurkace magnetic and®electric 
fields were recorded using the instrumentation as described 
below. The surface variations have been analysed to 
determine the local resistivity ‘structure of the "earthy >-In 
this thesis the polarization characteristics of the magnetic 


field variations will be analysed. 


Three mutually orthogonal induction coils were 
buried in the ground to measure the vertical (Hz), the 
geographic north-south (Hx) and the geographic east-west 
(Hy) components of the magnetic field variations. The coils 
Consist of 32,000 turns of #28GA°enamelled copper wireson a 
60 Mach xe cOr7s> MnehMten Mermeabt rity core eU-The -Fotal 
resistance of the coil is 845 ohms and the total inductance 


fs g50 Heriryetat ONL Hz! 


The: signals from the. colls was ted invova,pre- 
amplifier/post-amplifier system. The post-amplifier has 
gain settings from -20db to +20db in steps of 10db to 
permit the recording of various levels OL Signal, strengun. 
The signal was recorded on magnetic tape using @ seven 


channel Precision Instrument analogue FM tape recorder with 
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a tape speed of 15/16 inches per second. 


Six channels of data were recorded on the tape. 
Three channels record the mutually orthogonal magnetic 
iveldivarvautons two channels record the telluric field 
variations, and one channel records a WWVB time signal. 
The six signals were simultaneously recorded on a six 
ehannel Brush paper recorder for” visual editing prior to 
the: @ipitization:- process. A block diagram of the recording 
System is shown in Figure 2.1 and the frequency response 


of the magnetic system 1s shown in Figure 2.2. 


PAN Analogue to Digital Conversion 


The Brush paper records were examined to insure 
that they contained sufficiently. large amplitude signals, 
Momcisconvinuities , and did nov exceed the dynamic range 
Srmihesrecording apparatus (Saturation)... Suitable records, 
3.6 hours long, were digitized with a digitizing interval 
of 32/15 seconds and a consequent Nyquist frequency of 
0.234 Hz. A block diagram of the analogue to digital 
conversion system is shown in Figure 2.3. The analogue 
signal was passed through a low pass filter Wiuh a corner 
at.0.125 Hz and aliasing of the signal was thereby 


eliminated. 
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Bieure 2.1 


Block diagram of the recording system. 
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Figure 2.2 


Calibrated frequency response curve for 


the magnetic recording system. 
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Figure. 2.3 


Biogk diagram of the analogue to digital 


conversion system. 
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Digitized data were recorded on two, seven channel, 
asynchronous, incremental, Kennedy digital tape recorders. 
Bachydata point was storedtin asl2sbit binary form with 
Bece lenscisegniticant bitwuti lized ‘for ichnannel identifi-_ 
cation. Data blocks of 3600 two-character words were 
recorded alternately on the two Kennedy recorders to 
provide the Inter-Record-Gap required for input to the 
IBM 360/67 digital computer. The data were read from the 
two seven track tapes and stored on an IBM nine track tape 


for analysis purposes. 


EES epectral Analysis 


The autopower and crosspower spectral estimates 
were calculated in order to use the method outlined (see 
section 2.4) for determining the polarization characteris- 
Pece ames unnary -Or=tne~vechnigue used! or “obvaining*tie 
autopower and crosspower spectral estimates will now be 


given. 


easier Removal of Mean and Trend 


The spectral analysis technique requires that the 
time series originates from a stochastic process with zero 


mean value. lf the mean is nov zero, the power spectrum 
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Wit have a largce NC estimate which distorts estimates at 


Gunner sirequencles. Slr the record consists of N data’ points, 
Un» Vic Honest wa. eticmeolimave OL the mean is given 
by 
= 1 N 
ae N (bh eae 
ie n 
n=l 


This mean value was removed from the time series. 


Next, the method of Bendat and Piersol (1966) was 
employed to remove any slowly varying linear trends which 
Can cause distortion of estimates. These trends could 
Semoealifromeinsprument drift or actual trends. The average 
Slope of the time series u(t) can be defined to be, 


aL 
21/8 


a = zs Bi 


u T/73) (21/73 u(t)dt - ae u(t)dt] (252) 


where T is the total record length. Rewriting equation 
Peeing a Lorm culvepie Lor Cigivized, data) gives 
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where A is the digitization interval and v is the largest 
integer less than or equal to N/3. The data with zero 


mean value and zero average slope is then given 
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It is important to note that prewhitening the 
time series was not considered necessary since the micro- 
pulsation power spectrum exhibits a natural 1/f dependence 
(where f denotes frequency), which is counteracted by the 
output of the induction coils which is proportional to the 
frequency. Hence, the coils effectively prewhiten the 


Signal spectrum. 


eee Calculation of the Autopower and Crosspower 


spectral Estimates 


The autopower and crosspower spectral estimates 
were calculated using the method of Blackman and Tukey 
(1958) combined with the Fast Fourier Transform (FFT) 
algorithm of Cooley and Tukey (1965) as applied by 


Gentleman and Sande (1966). 


The FFT technique assumes that the data is cyclic; 
therefore, to avoid problems in the frequency domain, it is 


desirable that the initial and final data points of a time 
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Zeroes must be added to the time series to prevent 
Flew tGroguecea Cyclic structure ©f the data from distorting 
the autocovariance and crosscovariance estimates. The 
minimum number of zeroes required equals the maximum lag 
number GO be calculated. The FFT programs available require 
Liatevie TOvalenumber of points, N* 3 “penne analysed must be 
faclOlaylel arr BAe tediveowm 2s \herefore, the titel’ series 


or N' points 1s given by 
N' = N + number of zeroes (255) 


The discrete Fourier transform of the time series 


x is given. by 
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Weel Fen ee Nod (250) 


where Was the digital frequency =.ndex. | The angular 


frequency, w, is related to W by 


Sim larly, a. second vime series , ve will have a discrete 
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Fourier transform given by Y(W). The raw, unsmoothed 


crosspower estimate, Pxy(W), is then given by 
Pxy(W) = X*(W) © Y(W) C20) 


The asterisk denotes the complex conjugate. The 


Wiener-Khintchine theorem states that the crosscovariance 


and the. .crosspower spectrums are a Fourier transform pair. 


mherefore, the crosscovariance function, Cxy(t), is given 


by 


4 2TiWt 


Le aan. 


1 N-1 
Gx Cts T [ x {X*(W)-Y(W)} e 
This indirect method of calculating the crosscovariance 
requires much less—computer time than the direct 
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The Parzen window (Parzen, 1961) was used to smooth 


the raw autocovariance estimates. This window requires a 
negligible amount of computation time and gives positive 
definite spectral estimates with extremely low side-lobes. 


The Parzen lag function is given by 


R(t) = 1-6(4)? + 6(F)* |r] < 5 
Se po (=) 1% 3 <\t | <em (2.9) 
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where m is the maximum lag. 


The smoothed crosspower estimates were then 
evaluated by the following computation: 


non _2niwt 


Peet eCeaCT R(t ie, ites (2.10) 
tT=0 ; 





The smoothed autopower spectral estimates were obtained 


Piece Mularvfashionsby substituting 
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2o5e Removal*ofithewDigitizingyYPhase Shift 


Resmen lv bened tne sect ionic wes gsi x channels of data 
Mellencamp led sequentially Go give a digitizing interval 
Sumo 7 ioesecoids, Thus, the ldieitized points on two 
adjacent channels had a real time separation of 32/90 second. 
This iicroaquced a frequency dependent phase shittsin the 


crosspower estimates. 


Poscalculave tne acluad value, or thie phase echiti, 
consider the Fourier transform of a function of time, x(t), 
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> My larlyarsecond function, y(t), with a time translation 


Oligo. Jocethie tTollowing Fourier transtorm: 


Y(w) = ao 7 SHE Ne wat 
C2. 2) 
= et BY (yw), 
The crosspower is given by 
Pay Car Sere ee Per (ie PCP ey Ca) arene) (2.13) 


MremuGreaide y(t) wane Lect ced ands t xtt jlo. Samouled 
Picse wunenaG 15 zero and 8 represents the welavive, time 
difference between the two sampled channels. The cross- 
power spectral estimates obtained by the analysis of the 
horizontal orthogonal magnetic components were multiplied 
-iw8 


by e to obtain the correct crosspower phase estimates. 


Finally. 16 sheuld be noted that .inductiongcol 1s 
measure the time derivative of the magnetic flux $(t). The 
following property of the Fourier transform allows power 


estimates of the horizontal orthogonal magnetic intensity 
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For a coil, the flux is linearly related to the magnetic 
intensity. Therefore, when the calibrated instrument 
response was removed in this analysis, power estimates in 


Onitemot y¥*/cos Were obtained. 


Pal Determination of Polarization Characteristics 


The autopower and crosspower spectral estimates 
were used to calculate time varying polarization character- 
istics of geomagnetic micropulsat®fons in different ‘rrequency 
Baads -bynrtoLlowLngeithep»treatment+outlined by Fowler et al. 


(1967). 


Throughout the following, it is assumed that the 
vector disturbance recorded at a fixed point in space is 
quasi-monochromatic at theymean signal frequency, ft, Gi} 
interest. Quasi-monochromatie means] that 1f the irequency 
spectrum of a signal has a frequency peak of bandwidth AT 


centered at frequency f, then 














7 - teh 
j 1 
“y . _ 
:* ' 
ms  « ic : 
£F 
‘bayvet sd.ot YH bas xi 
oO WS \ ~~ P _twi- (t)o6 @ 
+5 s(t)b Lowe ti a 7 
x ~ ae 
(HL .S) 
= “~ 
(wo wt = 
oftengsm sit ot HDswsie7r ylusenil at xolt edt , {1 
tnsmuytent betsrdtiso ¢ aadw . sxyoterodiT v 
nf esvJemistes yswoq ,etaylens sid? at bBsvomet 2 
.benitatdo sxsw ego\*y to ad tay 
6 
eoivetistosisnd noftissitasf[ol te aobisatinreted - 
estemisves [attosqe iswoqezo%> bas teswododus sedaT 
~isJosisdo mottissiuslog sntvisy smist Siyolss ot beey 


YOMsypstT sasisitib afl enotiselyugoto im ofisnasmos3 to 


and a. Lr ‘r BR mow Ff heen ae r it ani * ‘ =: ‘ 
-[s 39 "si{wot yd Bbootlidso tremssit sdf anifwollot yd ab 


- 


Sit ted? bomyees ef ti ,antwoflot srg tnonguesdT 


uy 
ef sosqe of tatog Hsxit 5s ts hbsbtoos4 sons: disvsekb “i 


/. 
J 






an 






: to mj « YONSUPS TT fenate nsem sdt ts ottsmotdoonom=E. 
? 7 « 
7 —,* si¢ Tf jsd¢ ensom Sta BHee SOREN LaRtG taps 
7 


9A tb f: isin to Aseq wonoupart 8 esi [sngle s to me tt 


a a 7 : - mrs -_ y3 ) 









NEA oa js (25155) 


If Pxx(w) and Pyy(w) represent the autopower 
estimates of the north-south and east-west components 
Respectively Cf the horizontal magnetic field antlensity , 
and Pxy(w) represents the crosspower estimates at angular 
frequency, wo then 4 coherency M@vrux Can bevel ined: ae 

Bx Gon ey (iG ) 
JI(w) = (eee, LOD 
Pyx(w) Pyy(w) 
The matrix is Hermitian since Pxy(w) = Ee The matrix 
may be expressed as the sum of a coherency matrix, L(w), 
HoOeNate pOlerl zed portion ofp the signal and a coherency 


Matrix U(w) for the-unpoVarized portion: 
Jide = le Cue re CU ene) 


The polarized portion may be written as 


C(w) D(w) 
L(w) = 4 Con 3) 
D (w) E(w)|, 


and the unpolarized portion as 


F(w) G(w) 


U(w) = CenL9') 


et Aa 


For completely unpolarized signal, the mutual coherency 


must be zero: 
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le(w)| = 0, 


© (220) 
EXC pee aes Np Pale 


since the matrix is Hermitian: 
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if .ene component.were largerithan,the other, some degree 
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whe) coherency macrix,fornehe unpolarized portion of the 


Seenals can be rewritten: 
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Equation 2.16 can now be written as 


PxxCa)  Pxy Cw) LEI LOCKS) F(w) 0 
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Pyxtw) Pyy(w) D (w) E(w) 0 F(w) 

As» shown:in Born).and Wolf-(1959), the determinant 
of the.coherency+matrixefor,a.completely polarized signal 
ie always zero. This’ fact: enabiles, the value or I CAY) Saexe: 
be determined... Rewriting equation 2.22 yields 
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Taking the determinant of both sides gives 
F°(w) -— F(w)(Pxx(w) + Pyy(w)] + [J(w)| = 0, 


Where the vertical bars represent the value of the 


determinant. Therefore, the value of F(w) is given as 
1 

Bw eal Pex Ca) + Pyy(w)] + 4, (Pxx(w) + Pyy(u))* —- Ald(w) |). 
(2.24) 


Pees rOOn WIL the «pOS ttive sien is rejected=since 1t gives 
negative values! for’ the’ autopowers of the polarized, signal, 


CCm) and E(w), Bn equation! 2.123: 


Thesvedues Ofo¢(w)yeb(w), D(w) and F(w) can now be 
determined in terms of the original coherency matrix 


(equation,c. 16 )s 
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Bye Ucime sche "coherency Matrix fom the polarized isignal, 


Ee) the polarization ellipse may be deteymined from 
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Expanding equation 2.27 gives 
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Equation 2.28 is the equation of an ellipse in the coordinate 
system (x,y). The angle, 9(w), through which the coordinate 
axes must be rotated to align with the major and minor axes 
Simeone vellipse as called vine. ene levor polarization. _ This 
tosieo sc. by Shown ve be 
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Substitution of the values from equations G2.25)) allows 


calculation of 6 from spectral estimates. Thi see ves 
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If 6(w) is positive, the angle is measured clockwise from 
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Lvecomunowetiaipelrerty,. Lhe ellipticity 12s derined in terms 


of the angle 6(w) where 


1(D (w)-D(w)) 


sin 28(w) = ————— + .___\_ 
hCCH) =BGw) = oF AD <w)-Diw) 1-2 


where i =v -l . 


Substitution of the values from equation (2.25) 


yields 


ilPyx(w)-Pxy(w) ] 


sin 28(w) = 2 
[(Pxx(w)-Pyy(w))? + 4Pyx(w)Pxy(w) J? 


The ellipticity itself is given by tan B(w) which ranges 
frome estore circularly polarized Signal,’ to 0. tor linearly 
Volearucedvsigneal., lf Bw) is positive, the sense of 
polarization is counterclockwise when looking into the 


progating wave; the sense is clockwise if B is negative. 


PIS Implementation 


The Hx and Hy autopowers, the angle, the degree, 
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andes the sense, of polarizatden,fandizphe ellipticity were 
peebted ase tunctions of perlod»andyor local’ time.» The 

3.6 hours record was divided into sections of 818 seconds. 
Each section overlapped the previous section by 382 seconds. 
Overlapping sections gave smoother changes in the autopowers 
and polarization parameters as they were plotted as functions 
Orevulmes forma Derticular period. kKesults were! printed by 

the computer with local time’ as the abscissa and period as 
the ordinate. All the autopowers and polarization parameters 
were divided into intervals and assigned a value. This per- 
mitted contouring constant values of autopowers and polar- 


iMmapronmparameters as’ tunerions of local vime and pericd. 


Examinatiaen of the computer printout plots of angle, 
ellipticity and sense revealed small isolated "pockets" of 
eontucederesults.« These: tended tol eccur when, the degree of 
polarization Waslsmall.s lo Wacilitate contouring. whe values 
foraancle.vecllipticityo and) sense) were: lefts blank»on) the plots 
when the degree of polarization was less than 50p. “sinces the 
signal was predominantly greater than 50% polarizeds this 
resulted in only isolated pockets of blanks. When contouring 
was being done, these blank pockets were ignored for ease of 


interpretation. 
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CHAPTER III 


THE RESULTS OF ANALYSIS OF MICROPULSATION DATA 


shad Recording Stations 


During the summer of 1968 measurements of the 
surface magnetic field a ee hes made at the 
University,of Alberta Geophysical Observatory* near Leduc, 
Miperce sone 5S 2 3 Nk = isos Wie eCOnLinuous recording 
was carried out between May 20 to June 20, 1968, and 
between August 1 to September 13, 1968. No recordings 
were taken between 0900 to 1000 hours LT** and between 
2100 hours to 2200 hours LT** as this was the time that 
tapes were being changed on the recorders and calibrations 
of the instruments were being made. Records that contained 
Mo Activity or that had: exceeded the dynamic. range of the 


tape recorders were discarded. 


Rankin and Reddy (1969) have stated that 
anisotropic. electrical conductivities exist Jmepne Pre 


cambrian basement in the vicinity of the Observatory. It 


* For brevity, the University of Alberta Geophysical 
Observatory will be called the Observatory. in: the 
remainder of this thesis. 


*% Local Time (LT) is Mountain Standard Time which as 
seven hours earlier than Universal- Time. 
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was felt necessary, therefore, to compare the polarization 
parameters obtained at the Observatory to those made 
simultaneously at other sites. Thirteen field locations 
were chosen in, central western Alberta and in southern 
Aigcela eee neLlr Locations. ere shown in Picure, 3.1. “Also, 
their geographic coordinates are given in Table III. The 
dates and beginning times of the records analysed are given 
in Table IV. The records last 3.6 hours. Similar equip- 
ment as that used at the Observatory was mounted in a 
PeUCKe LO. recordsay tne, fieldveitess When two or three 
Suitable eleven hour records had been recorded, the truck 


was moved to the next site. 


Locations were chosen on the basis of isolation 
frommenunan noise Contamination. Thus sites were as. far as 
possible.from: busy highways, electrical sources, construct— 
fon and oll wells. Alsosthe surface. relive? was’ small for at 
feast. several miles in-all directions from the eee aeaae 
station. These conditions were easily met on the Prairies 
of southern. Alberta. However, the stations of Nordegg, 
Foothills and Gregg Lake were located in the foothills of 
the Rocky Mountain Range and the surrounding relief was 
considerable at distances greater than half a mile from the 


site. 
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Figure 341 


Map of Alberta showing the. location 


oiP the recording stations. 
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TRB Uap el 
GEOGRAPHIC COORDINATES OF RECORDING SITES 


STATION LOCATION. GEOGRAPHIC GEOGRAPHIC 
LATITUDE LONGITUDE 


Ader Flats Fee (tN LTA 57 14 
Camp Creek BHO 3'N TTA 32 
Carmangay 50° 5" N 1ay39 ke Wi 
Carrot Creek bce 6. N 15°50) 4 
Fioo Fidelis 53° 0e UN 116°48 'W 
Fox Creek 54°26'N 116255 'W 
Gregg Lake 53-307 N Vesa OG 
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Horse Thief Canyon 51° 34 Nn TLS 55 'W 
Nordegg Boo eo LN 116°05'W 
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ae Record Characteristics 


Since induction coils were used to record the 
variations in the surface: magnetic field, it-was difficult 
pome lass iry DpuUlSeulons aato Pits sand re’ s. dhe ampiitudes 
of the various period bands were out of proportion since 
dH/dt was being recorded. However, a close examination of 
the records indicated that the pulsation events between 
2200 hours to 0500 hours LT could be classified mainly as 
Pi events and those between 0500 hours to 2200 hours LT 
were mainly Pc events. Examples of the analogue records 
for the north-south (Hx) and east-west (Hy) components of 
the magnetic fleld variations; are shown fae eUuresms .< 


EVa0) Geen 


The vertical component (Hz) of the magnetic field 
variations was small when compared to the horizon an 
components for almost all the records. Any appreciable 
activity usually was confined to periods greater than 300 
seconds. Therefore, the polarization ellipses were assumed 


to lie mainly in the horizontal plane. 


The Kp indices are also given ia Tables iy. selec p 
index is a 3 hour index which is designed to measure the 


"planetary" variations in geomagnetic activity. It has 28 















eotieliedosiend breve , te 

sit bicsss o¢ bees stow alton nottoubak sonle 
$ivoeltitb saw dt . bi ott oftetissm sostuve antt ne snot¢sitsv 
aebudtiqmea oAT .2'oF Bas ett ogmt enotiselug ytteests ot 
sonte noitiogo1q to tuo sisw abasd botreg evotuey ort to 
to noftsninsxs seolo e .asvewoH .hsbaosss sited esw JO\Hb 
xsswied etnevs nolttselud edd tsdZ betsobpnt abtooss edt 
28 are beltlees{o sd bluos Td awed) 0020 ot “etuond O0SS 
Ti siyon OOSS oc etuod 0020 nsewted seodt Bas siaeve 2B 
abioosy sugolsns edt to aslamexe .edtriave of yintam stew 
to etnenoqmos (yH) teew-tase Bas (xH) fitjvos-dAtton sd¢ 10% 
S.€ astugit§ af mwode ors enottskasy bLStt oftengsm ed? 
| £.€ baw 

Hfett oftjengsm sdt to (sH) tnsnogqmos Sas isiey sciT 
fsinositod sdt o¢ bsisqmos oorlw' fieme asw sii: sted 
eldsisesigqgs yok .ebrooss sat {{s teomis to? a 
QO€ asdt tetse1rg shottaq o¢ Genie tase cial yilsveu yivisos 
bemsees stew aseqtii[s noftsstrslog sit Sto TStoAT . abtioose 
-onsg {stnositor eat mt al ett of 


ME — ot revs cals os & 


7 - “LL. 35 et it Ww Sa 


J Ue ’ : », os t 
_ Pas ™ Saf 7 Si. wT. 
ck J 28 ai 


a = 





7; 


Pieure: Sez 


Analogue records for May 30, 1968 as 
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Analogue records for August 17, 1968 as 


recorded at the Observatory and Carrot. Creek. 
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grades ranging from 0,., which denotes exceptional 
quietness, to 9, which denotes the most intense storm. 
As can be seen from Table IV most of the records used 
were made at times when the Kp index was 2 or greater. 
Two Kp indices are given if the record spanned two Kp 
three hour, index values. Thus most.of thesrecords used 


represent times of reasonably disturbed magnetic conditions. 


4.3 Comparisons Between Simultaneous Field and 


Observatory Results 
a. 30 Individual Records 


A careful comparison between Observatory and 
simultaneously recorded field station power spectrums 
showed that they were similar in all cases. Examples of 
field and Observatory power spectrums for the north-south 
component (Hx) are shown in Figures 3.4 and 3.5. These 
power spectrums have not had the instrument response 
removed. Therefore, they are essentially the prewhitened 
spectrums and the contours shown are used only to deter- 


mine the periods at which there were power peaks. 


As mentioned in Chapter II, plots of polarization 
angles were made as functions of time of day and of period. 


Examples for the field and Observatory are shown in Figures 
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Figure 3.4 Prewhitened digital sonograms for May 30, 1968 


for @Ghe pObservatcory and Horse Thief Canyon. 
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Prewhitened digital sonograms for August 17, 


1968 for, themObservatory and Carrot Creek. 
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3.0 and 3.7. No consistant visual correlation could be 
found between these plots and their corresponding power 
Spectrum*@plots. That is. the polarization angele did not 
seem to depend on the power present. The angle could 
remain constant with respect to time whether there was a 
power spectrum peak present or not. On the other hand, 
the angle could change rapidly irrespective ofthe form 

OF theppower spectrum plots “This could Suggest’ that: the 
polarization parameters are dependent on local ionospheric 
Genditiens above the recording station and, not on,the 
source mechanism in. the magnerosphere... The source 
mechanism, then, would be responsible for the amplitude of 


The pulsation event. 


OUbeot*tnemw4>o Dalts orm. eldsandsOpservatery 
records analyzed ..j pairs, snowed @digvle ortno redationship 
Go, each Other instheir angie of polarizations plovs. 
However, angle plots recorded at the same field station 
at.a different time would show a cocdrscorrelation with the 
Opservatory This" indicates! that themiack olf commelarion 
is not entirely due to geology of the field station. If 
one remembers that the power spectrum plots were found to 
be similar, then this result.is probably an indication that 


while the macrostructure of the pulsation event remains 
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Figure 3.6 


Angverorpo@tarizataon—for-May-30'- 1968 


for the Observatory and Horse. Thief Canyon. 
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Figure 3.7 


Angele of polarization for August 17, 


Tor the Observatory and Carrot Creek, 


1968 


Da 


ad 





=. 


i 


Seer . VL teuguA rot noltssiysiogq to sigad YT itt 


-Ass710 tor1s0 bas ytotsvisedO add rot 


oN Bary 


"3 
>. (eur 





ANGLE OF POLARIZATION 


AUGUST 17, 1968 


\\ 
LS WVU 


STL TUE 
— \ 


CS 2 # 








(SGNO93S) GOlu¥ad 


LT (HOURS) 
CARROT CREEK 


LT (HOURS) 
OBSERVATORY 


ean | 
‘MOIT: ASIAAJOT 


- “Ba agt et 31 Tay 





bE #IOD. (2@C0noe) 





; -? » 
4 * : 
~—e a 
7 
iz | 
: ~ 
h, _ 
j ) 
i = 
| ene. 
» 
~ 
f 
t > 
_ =. 
- 
| j ‘ nw 2 
: : . * = 
-] 7 = - = gape AL od 
a — 
4 - n 
E ; i? 
—— ™ 
rr } f 
t i bi Fi TJ 
x ; . 
les 
o 


54 


relatively constant over a distance of several hundred 
Kplometers, the phases’ of the signal can change. This may 


Dossibly be due to Tonospneric inhomogeneities. 


Somewhat the same results appeared when the sense 
of polarization was considered. Figures 3.8 and 3.9 show 
examples of the sense as determined from field and 
Observatory records. Five field and Observatory pairs 
showed little correlation. That is, the sense would switch 
from clockwise to counterclockwise polarization indepen- 
denclysay eacn station: Once again, the sense of polar— 
igation plots showed little, if any, relationship to, the 


power spectrum.plots. 


St Sipe Averaged Results 


While each angle of polarization and sense of 
Lomarizar lon plot nadhvery interesting characteristics, 
nothing definite could be said about these polarization 
parameters in general. Therefore, all the plots were 


averaged to determine more general characteristics. 


Ca) Sense of Polarization 


The sense of polarization obtained from the 


Observatory records were averaged to determine the 
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Figure 3.8 


Sense of polarization for May 30, 1968 


for the Observatory and Horse Thief Canyon. 
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Figure, 3.9 


Sense of polarization fortatrust 17, 1968 


for *the™Observavory and Carrot. Creek. 
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dominant sense as a function of period and time of day. 

All the plots were sorted into hour intervals according, £o 
the time of day at which the records were taken. This was 
done for each period estimate. (As was shown in Chapter II, 
the Fast Fourier Transform technique gives discrete auto- 
power and crosspower spectral estimates at certain periods. 
The polarization parameters were calculated for each of 
these spectral estimates. The term, period estimate, means 
the sense of polarization that was calculated from the 
spectral estimates.) Then, the number of times that the 
polarization sense was counterclockwise, was divided by the 
total number of polarization sense. estimates to give the 
percentage of sense estimates that were counterclockwise in 
an hour interval at a given period. Figure 3.10 shows the 
result which was obtained. The total number of sense 
estimates which were used in each hour interval are dis- 
played at the bottom of Figure 3.10. The segmented sections 
of the contour lines represent the times at which no data 


were available. 


The figure shows that all micropulsation events in 
the period range from 80 to 273 seconds are 60 to 70% 
Dolarized in the counterclockwise sense, except between 


1300 hours to 1800 hours LT. The intrusion of counter- 
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Figure 3.10 Percentage of counterclockwise polarization as 
NR DUneD 1 enwToT—Timeror oayoand—pf period.| Alt 
of the Observatory records were averaged. to 
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clockwise polarization into shorter period events between 
0500 to 1200 hours LT supports Zybin (1967) who found that 
the sense was 80 to 90% counterclockwise for all micro- 
pulsations in the period range from 30 to 1200 seconds in 
the local morning hours. In the period range above 100 
seconds), Pigureo3s. 98showed that 80 to 90% of the micro- 
pulsations events were polarized in the counterclockwise 
Selioce Neoie Oca Lomornings bours 1. This Wailso Vs5 consistant 
with the results obtained by Rostoker (1967b) who found 
that 70% of the Pi2 events had a counterclockwise sense 

Of polarization. If one remembers the predicted inverted 
sense between the northern and southern hemispheres on 

the basis of Alvén wave propagation, then, this also 
agrees with the results of Christoffel and Linford (1966) 
who found that 61% of the Pi2 events recorded in New 


Zealand had a clockwise sense of polarization. 


Frequencies which correspond to the Pil band have 
LOmvOs SO ce On their contributions. polarized. in, the clockwise 
Girection_in,thbe Jocal,nigshirime hours, This clockwisessense 
@onvinues., untad. 1100 or 1200 shoursmaLl and, then changes 
abruptly. They then become 60 to 80% counterclockwise polar- 
ized. Campbell (1967) states that during the daytime max- 


imum occurrence period of Pc3's, the counterclockwise sense of 
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polarization predominates. The maximum occurrence period 

has been found by Jacobs and Sinno (1960b) to peak at 1200 
Po. ¥1830 Oftouns CRN nlenertiern Imirddke bavittudesito The ire savbts 
shown in Figure 3.10, then, agree well with Campbell (1967). 
The sense of polarization corresponding to the Pc4 band is 


Seo me hyemiomeréc oni wsied 3 


Campoell alse staves vat, tre angie of 
polarization undergoes a diurnal change with which the 
sense also changes. No conclusive evidence for this was 


homndmacewill pe seen slim Dart .b) "Oreunis section, 


Cb) Mnple of Polarization 


Examinattron of individual plovs of angle of 
POlarizacvion versus period and time of day, showed that it 
was difficult to determine general trends. Therefore a 
computor program was written to compute the average angles. 
As was done for the sense of polarization, all the 
polarization angle estimates from the Observatory records 
were sorted into hour intervals according to the time OF 
day atv which they were recorded. Then, the angle estimates 
were grouped into 5 period bands which were: 10 to 20 
seconds, 20 to 45 seconds, WeetG G8 ceconds, oc toeis0 


seconds, and 136 to 273 seconds. The angle estimates for 
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each period were then averaged to give the average angle in 
a given period band. These averaged angle estimates over a 
given period band were then further sorted into angle groups. 
There were 12 angle groups ranging from -90° (i.e. geographic 
Wesc)yctomtToO0 (tee. geographic, east}! in increments of 15°. 
Maem pemcentace or themangiles Talisima in each otethesl2 
eneleseroups Was displayed on histograms and are shown in 
(ape omom Le onrougn Figure 3.15. sFrieures 3.11a andes. 11b 

are for the period band 136 to 273 Saco. 7 Lae andr 3% l2b 
Pom iocmco 136 seconds, 3.132 and 3.30 for 45 to 60 seconds, 
Betis cand 3.140 for=c0r vow > seconds, .and, Figures 3-1 5eanand 
25955 areidorgtbhespertodhbandg1® To 2Oisecondson Therrecords 
from the 13 field locations were averaged together in a 
Similar manner and the results are displayed beside the 
Observatory histograms) for comparison. It must be 
remembered that each field station record corresponds to 
armeOovsenvavory record. ~iheretores any differences between 
the field and Observatory histograms are mainly due to 
changes in geoelectric structure of each field station and 
changes in the signal as it propagates from the source to 


the two recording stations. 


An examination of Figure 3.11 through to Figure 


Seo, Lom che Observatory, immediately reveals one interest- 
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P#@euressilla and Figures sie 


Histograms for the angle of polarization averaged 
over hour-—intervals for the Observatory and field 
records. These figures are for the 136 to 273 


second period band. 
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Fisuress..ca and Figweems.12b 
miscograms ror the angle Of polarization averaced 
over hour intervals for the Observatory and field 
records. These figures are for the 68 to 136 


second period band. 
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Figures .13a and Figures3. 13b 
His Lopramsetor. tne "angle “otepolarization averaged 
over hour intervals for the Observatory and field 
records. These figures are for the 45 to 68 


second period band. 
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Figure 3.1l4a and Figure 3.14b 
Histograms for the angle of polarization averaged 
over hour intervals for the Observatory and field 
records. These figures are for the 20 to 45 


second period band. 
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Feyeure=3el5a and Figure 3. 15b 
Histograms—ror=chnesangle fofepolarization averaged 
over hour intervals for the Observatory and field 
records. These figures are for the 10 to 20 


second period band. 


66 


we Daa 





aoe 


“st 1s (AEWA 


» : 
défl.€ otuzati bas scL.£ 9 
begsisve notisstusfog to efgns eit tot emstgovetH 


“Sa 
blett bas yrovsvreadO st rot elgviednt wod eve. a 
OS oF OL oft yO. ots eougtt seediT -abtoss1 | = 
.basd bolseq Bhooge” = 
oo 
ar G eS 
a) 
” zt = 





IN EACH 15° INTERVAL 


arian 


PERCENTAGE OF SAMPLES 


OBSERVATORY | FIELD 
10-20 SECONDS 


yo | 9000 yo | 0000 


0100 


0200 


0300 


0400 





0500 


0600 











: 0700 | uo 0700 
20 20 
0 0 
<4 0800 - yo 0800 
20 20 
0 0 
a4 1000 . yO 1000 
20 20 
0 0 
TT) 1100 40 1100 
:: al: : si 20 
0 0 
-90 -45 O 45 90 -90 -45 O 45 9 


0 
ANGLE OF POLARIZATION (DEGREES) 








| 
7) 


- 


; iv Ne ic a 
. ne 


gj314 ¥AOT “AVA 280 












: gvuoos2 OS-Ol 
‘ff ¢ = 7] : us : . ; 
OF nad ee -. 
7 4 ane, _ 7 X ay ae 
. “a y : 7 J F . Q at a = ™ m9 pr ate ; tT é 7 
OC or 0010 bat i ,) 7 
>| ’™\ * v2) ‘ i 9 . : 
" - - ; 
[ ~~ v A. ¥ 
f { , 
; 7 
ic ’ 
: 5 - 
; ; + fe bes | 0 _— ‘ 
, 5 - 
a“ im ° e ’ ; 7 0080 1) tl 
| | os ta 
a ‘oS t , 
TT : . ’ rt . a we ~ ere - 
~ | ——— a Co 
; 0020 we 
Os 
Ne cs 
.. a J —— oa 0 
| —_ 
| 180 F ow 
4 
wv zi oS 
. eaall 4 on; § 4 . 
’ ’ = b . 
} a : 
J 
ws Eo ; 
an 
— os 
q f 7 >; be oe a) 
us ' - RN Ee 


0080 


TT 


| 
' 






HTT 


PERCENTAGE OF SAMPLES IN EACH 15° INTERVAL 


OBSERVATORY 


lO-20 SECONDS 


= 1200 
20 
0 
uO 1300 
20 
0 
= 1400 | 
20 
0 
= 1500 
20 
0 
4 1600 
20 
0 
uO 1700 
4 Hier: 
0 
a4 1800 . 
20 
0 
ug 1900 | 
20 
0 
uQ + 2000 
20 | 
0 
Te) 2200 
= mi ee 
0 
uo-| 2300 
20 
0 
-90 -45 O 45 90 














FIELD 

up 4 1200 
20 
0 
Bah gas 
20 
0 
4“O- 1400 
20 
0 
spe 
20 
0 
ae 1600 
20 . 
0 
up J 700 
eer 
0 
Yo 1800 
20 
0 
uo 1900 
20 
0 
up 4 2000 
20 
0 
up 42200 
ee 
0 
SE Te upeens 
20 
0 

-90 -45 O 45 90 


ANGLE OF POLARIZATION (DEGREES) 








ak 2 one of 





67 


Migstecu, Ne Mayoreaxis Of Une polarization eltipse, falis 
predominately between 0° to 60° east of north. This result 
COulG be explained am tnree Ways. “lt could be*due* to 
geological structures situated in the near vicinity ofthe 
Observatory. However, an examination of the corresponding 
field histograms shows that they also have this general 
trend. Therefore, “1ibthezaneclie oftpelbarnhzatbionals bheineg 
influenced by geologic structure it must be a much more 
extended feature.» The Rocky Mountain Range makes an: angle 
of roughly 45°W of geographic north (see Figure 3.1). 

Zybin (1967) made the observation that the major axis of 
mae polarizavion ellipse tends to be perpendicular. to the 
Strike of. major geological features. This might explain, 
then, why the field and Observatory histograms display an 
angle of polarization lying between 0° to 60°. The third 
explanation 16 that the polarization angle is aligning with 
the magnetic meridian. At the Observatory, the magnetic 
Meridian ts 23° E of geographic north. Several workers 
(i.e. Campbell (1967)) have noted that polarization ellipses 


have vended to be aligned in this fashion for Pce5 events. 


A closer look aL the histograms revealed more 
features p),n1) Of therhisvograms Lor the Observauory svucnced 


TO Wlaveran angle of polarization in the 15° to 60° range in 
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thes loon bani ehteroums, G2s00shours to 0500 hours LT). The 
histograms seem to become more spread out and scattered in 
the hours between 0600 hours and 1100 hours LT. This 
corresponds roughly to the. intrusion of counterclockwise 
polarization into lower periods as is displayed in Figure 


Gre le 


Dn Rieures 3 wlkbseGe ei ct O-<l 36 eeaond period: band). and 
Figure 3.12b (136 to 68 second period band) the peaks in 
the histogram seem to lie in the 0° to 30° range between 
1200 hours to 2000 hours LT. Between 2000 hours and 2300 
heuneamian eUney swine to their nighttime position of 15° to 
OOee, 

The histograms representing the 68 to 10 second 
period pulsations in the time interval between 1200 hours 
and..2000,.nours, LT. ane, more. comp lexi. This complex form is 


ale o. evident.in.Figcune. 3.10) for, the..sense of polarization, 


Thus, it seems that for the longer periods (68 
seconds to 273 seconds) the angle of polarization lies 
between 15° to 60° in the local night hours. In the 
morning hours (0600 hours to 1100 hours LT) the situation 
is confused. However, between 1200 hours and 2000 hours LT 
the angle is between 0° to 30°. Then between 2000 hours 


and 2300 hours LT there is a clockwise swing in the angle 
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backoto 1ts nighttime diréciiom 


(c) Degree of Polarization and Ellipticity 


In Chapter II the degree of polarization was 
defined as the ratio of the polarized signal intensity to 
the total signal intensity. An averaging process, similar 
coOmene ONC IUSed EO OOtTaln Pigure 32.10 was carried out. The 
result showed that on the average, the signal is 65 to 75% 
Polar 7eGmnorwal peri ods, and at all times of the day. 
Similar results were found from an analysis of the field 


records. 


Again, the same averaging process was applied to 
Ehesellipticity vazues obtained from Observatory: records. 
It was found that, on the average, micropulsation events 
Hevesaurecvio. of minor to major axis of 0.2 to.0.4.- Thus 
iG can be expected that polarization ellipses will be quite 


elongarcds) similar results were found for s.the field) dara. 
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CHAPTER IV 


CONCLUSIONS AND SUGGESTIONS FOR FURTHER RESEARCH 


Mrs Conclusions 


mir unis work, methods were déveloped to analyse tne 
polarization characteristics of geomagnetic micropulsations 
Ugrlizing Taree emouni.s,Of data. Whittle results from*@indivi—- 
dual records were superficially interesting, more conclu- 
Sive trends were obtained by using statistical averages. 
The individual records did show however, that very similar 
polarization results could be obtained’ from Simultaneous 


field and}Observaboryurecords: 


The averaged sense of polarization values indicate 
(hat Micropulsations in the period*range* from’ 30" to 273 
seconds were predominantly polarized in the counterclock- 
wise sense except between 1300 hours to 1800 hours LT. 
For micropulsations with periods below 30 seconds the 
polarization sense is different than for the longer periods 
in that there is a semi-diurnal variation in the sense. 
The sense goes from being predominantly clockwise in the 
local morning hours centered at 03:30 to strongly counter- 


clockwise in the afternoon hours centered at ) Bo 3.x Lt 
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may be further noted that these times correspond to those 
in*which the longer period micropulsations are-less 


strongly counterclockwise polarized. 


The angles of polarization were also averaged and 
the results showed that the bulk of the angles lay between 
0° to 60° Gast. of north. The strong similarity between the 
field and Observatory stations indicates that local geology 
De tnotnascominantohactorlinadeterminingoatheranghesof 
polarization. Thus the observed angle of polarization may 
be attributed to the source of micropulsations or to 


regional geologic trends (such as the Rocky Mountains). 


The average angle in the local night hours lay 
between 152 towseosseastlbefonorthvallinithe/bimetbetween 
1200 hours to 2000 hours LT the angle lay between 0° to 
302% 

The records used in this analysis were not 
preselected except for a minimum level of acrivatyyvas 
described in Chapter III. While nighttime records were 
considered to be predominantly of the Pi type and daytime 
records predominantly Pc, the results shown here tadicaue 
that no significant difference exists between these types 
of micropulsations as far as the parameters measured in 


this work are concerned. Indeed, spectral analysis could 





eaord 0+ bnoqasvron esmid esedd Sst beton asdtiw? am 
eeel sis anotiselugotoim botyeq reyn0l. odd. dotdw: at 4 
.bestislog ealwaooloretayos ulsnotte 


Sas beystevs coals o1tew nottestteloq to eslans eAT fat ¢ 


| 
; 
| >t | 
msewted ysf eolgns sit to Alud sd? ted¢ bewode ativast edd . 
sid neswied yitisitmie gnoite efT .diton to dass 08, oF, °0 | 
wolosa {no0l todd astsotbat enotsste yrodevised0 baa belt 
to efsie edd antnimrsted at sotos? Jasnimob s tom at | 

ysm coléssitalog Io elgns bevisedo ed? apAT nolsssiuslog 

oF tO anotisalugotoim to somos sdt of betudixtss. od 


-CemistnwoM yaooh edd 2s dove) ebnext otgcolosg Lenobget 


ysl emod tigin [soot edd mt efgns egsreve Ty. ees 
meswied emit oft mI .dtaom to tess °93 oF LE noowded 
of “0 neewied ysl oelgnse oft Tl ewod 000S ot axyod ooss 


‘og eae 


ton stow ateylsns etdt at boey ebacoer Te be vate 
as UWwivitos to fevel mumtnim s t0% tqeoxs percolated 
Stew ebtooe1 ombitigtn eltdW IIL «sd¢qedo at bedttogen. fe 
cle eiantipa $8 eft Io yLtnsntmoberq ed alla | 
938 fe: A wore siueey ext OT Ultnsn ton 
ne Som a eageERstte 


——_ 


» ? " 


week ay = airplanes 


x 


Pees on 





te 


not be expected to resolve these categories. Therefore, 
the averaged results presented in Chapter III are really 
the. overalitpolarizationecharacteristics of themmicro= 


pulsation spectrum throughout the 24 hour day. 


Nee Suggestions for Further Research 


Since polarization characteristics correlate well 
between simultaneous field and Observatory records, it 
would besuseful to compare the same. characteristics 
obtained from micropulsations recorded at.conjugate,.point 
stations. Averaged results would perhaps clarify the mode 


of propagation of signal between the two hemispheres. 


The.data as it becomes available.from the 
Observatory can be used to establish seasonal and yearly 


TRends sinschespolarizationycharacteristicss 


Record sections longer than those used in this work 
willealso allow. the extension of this type.of analysis to 


longennperirods.: 
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